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Abstract. We report on the combined X-ray and ra- 
dio observations of the type Ic SN 2002ap, using XMM- 
Newton ToO observation of M74 and the Giant Metrewave 
Radio Telescope (GMRT). We account for the presence of 
a nearby source in the pre-supernova Chandra field of view 
in our measurements of the X-ray flux (0.3 - 10 KeV) 5.2 
days after the explosion. The X-ray spectrum is well fitted 
by a power law spectrum with photon index a — 2.6. Our 
results suggest that the prompt X-ray emission originates 
from inverse Compton scattering of photospheric thermal 
emission by energetic electrons. Radio observations with 
the GMRT at 610 MHz (8 days after the explosion) and 
1420 MHz (70 days after the explosion) are combined with 
the high frequency VLA observations of SN 2002ap re- 
ported by Berger et al., 2002, and the early radiospheric 
properties of SN 2002ap are compared with similar data 
from two other supernovae. Finally, the GMRT radio map 
reveals four other X-ray sources in the field of view of M74 
with radio counterparts. 



Key words: Supernovae:individual:SN 2002ap; Super- 
novae: general; Radiation mechanisms: non-thermal; cir- 
cumstellar matter 



1. Introduction 

Supernovae are known to be explosions of massive and 
intermediate mass stars, although the nature of the pro- 
genitor star for different supernova types, remains an area 
of long-standing research. Currently, it is believed that 
type II and type lb or Ic supernovae arise from core col- 
lapse of massive stars, while the more homogeneous type 
la supernovae are the results of thermonuclear explosions. 
There is a considerable variety among the core collapse 
supernovae spectroscopically and in their kinetic energy. 
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The 'hypernova' class among these are believed to have 
explosion energy significantly in excess of 10 51 erg, found 
in "normal" supernovae. The spatial and temporal near- 
coincidence of type Ic SN 1998bw with the long duration 
GRB980425 has provided further impetus to observational 
and theo retical studies of type Ib/Ic SNe and their pr o- 
genitors ( MacFadyen et al., 2001 , Meszaros et al., 1993|) . 

During a supernova explosion the interaction of the 
outer parts of the stellar ejecta with the circumstcllar mat- 
ter gives rise to a high energy density shell. X-ray emission 
is expected from both shocked circumst cllar matter and 
the shocked supernova mat ter (see e.g. |Chevalicr, 1982 



jClievalicr fc Fransson, 2001 ). In addition, the interaction 



region may also accelerate electrons to relativistic energies 
and amplify pre-existing magnetic fields which gives rise 
to nonthermal synchrotron emission seen in many super- 
novae. The radio and the X-ray emission give information 
on a region of the supernova which may be far removed 
from the optical photosphere (which has a smaller radius), 
although the conditions in the optical photosphere may 
determine the X-ray emission characteristics in some in- 
stances. In general, X-ray and radio observations in early 
stages of a supernova can be used to determine (1) the 
total mass lost from the pre-supernova star before explo- 
sion and (2) constrain various physical processes leading 
to X-ray and radio emission. 

The type Ic SN 2002ap was discovered on Jan. 



29.4, 2002 (Y. Hirose as reported by |Nakano et al., 2002 
) in NGC 628 (M74), at a distance of only 7.3 
Mpc. Based on spectral analysis of the early ob- 
servations, the epoch of explosion was estimated at 
Jan. 28.0 ± 0.5 2002UT QMazzali et al., 2002| ). For 
the purpose of the analysis presented in this pa- 
per, we will regard Jan. 28.0 as the date of explo- 



sion. The broad spectral features ( 


Kinugasa et al., 2002 


Mcikle et al., 2002, 


Gal- Yam et al., 2002), and a subse- 



qucnt modeling of its spectroscopic and photometric data 
( Mazzali et al., 2002 ), suggested that this was an ener- 

10 51 



getic event with explosion energy E 



10 x 
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erg. In this paper, we discuss our analysis of the XMM- 
Newton observation in the X-ray (O.f-f 5 keV) bands (Sect. 
||), accounting for the presence of a nearby object in the 
pre-supernova X-ray field, observed earlier by the Chan- 



dra X-ray Observatory (Sect. 2.1). We also observed SN 
2002ap in the 0.61 GHz and 1.4 GHz radio bands (see 
Sect. ||), and the implications of the GMRT upper lim- 
its in the context of VLA observations at the same epoch 
(Bcrgcr ct al., 2002) are discussed in Sect. ||. 



We also summarize the explosion parameters that we 
derived from the optical observations and modeling re- 
ported by Mazzali et al., 2002, as this is used as input 
for later sections (Sect. |4|). We have also combined the 



GMRT data with the VLA data ( [Berger et al., 2002| ) to 
derive conditions near the radiosphere. A combined analy- 
sis of the early X-ray and radio observations is presented in 
Sect. |5|, which attempts to constrain the multiple physical 
processes (thermal and non-thermal) that are responsible 
for the early X-ray emission. We also compare SN 2002ap 
with 2 other SNe (1998bw (Ic) and 1993J (lib) for which 
early multi frequency data is available. Finally, in Sect. ^| 
we discuss these results in the context of the presupernova 
star and its evolution. 



2. X-ray observations and data analysis. 

SN 2002ap and its host galaxy M74 (NGC 628) were ob- 
served by XMM-Newton, fo ur d ays after the discove ry of 
the SN. [Pascual ct al., 2002| and |Soria fc Kong, 2002] have 



reported the EPIC-PN detection of the supernova, as well 
as several other sources in the host galaxy. Here we present 
a spectral analysis of the XMM-Newton EPIC-PN and - 
MOS CCD detection, taking into account contributions 
from other pre-existing sources present in the extraction 
circle, near the supernova. 

XMM-Newton observed the field of view of SN 2002ap 
in the full-frame, thin filter mode, for the EPIC-PN and 
the two EPIC-MOS cameras. Simultaneously, it also ob- 
served M74 with the Optical Monitor in the UVW1. The 
full exposure in the EPIC CCDs was ~ 34 ks between Feb 
3.00 to 3.40, 2002 UT. Since Chandra X-ray observatory 
also observed the same field on June 19, 2001 and October 
19, 2001 for a total of 47 ks, we present our analysis, based 
on both the presupernova and post-supernova exposures. 

Early analysis of the XMM-ToO observations of SN 
2002ap showed that the supernova is rather faint in X- 
rays with a flux of ~ 10 -14 erg cm~ 2 s _1 . Because the 
psf of the XMM EPIC-CCDs is rather large (~ 50" to 
enclose 90% of the energy at 10 keV), it was necessary 
to verify the absence of, or to account for the presence 
of any other sources, however faint, within the spectral 
extraction region on the EPIC CCDs. 



2.1. Chandra observations of the presupernova field: 

In order to check if there were any other sources present 
in the vicinity of the optical position of SN 2002ap (RA 
= 01 ?l 36 m 23^85, Dec = +15°45'13'/2), we examined the 
archival Chandra observations of M74, carried out on 19th 
Oct. 2001. In the net 46.2 ksec of exposure time with 
Chandra- ACIS, we did not see any source at the position 
of the supernova. However, we did find a bright source at 
RA 01 ft 36 m 23!404 and DEC +15°44'59'/89, which, being 
only 14.9 " away from the SNe, was well within the region 
of XMM spectral extraction. Comparing the positions of 
bright, point-like sources seen in both the Chandra and 
XMM fields of view, we found that the relative astrometric 
shift between Chandra and XMM is at most ~ 5" . Thus, 
this object (hereafter, CXU J013623.4+154459) would 
have been detected within the 40" circle used to extract 
the XMM-EPIC spectra. 

Level-2 pipeline processed Chandra data was used for 
spectral analysis, since a check on the calibration files used 
in the archived level-2 processed data showed that no im- 
provement would be made by re-calibrating the data. The 
spectrum was extracted using the CIAO software, and the 
spectral analysis was carried out using the SHERPA soft- 
ware. The source was located on ACIS-S6 chip, and a 
light curve was extracted from a source free region of the 
ACIS-S6 chip to check that the observation was not con- 
taminated by background flares. We also extracted a light 
curve for the source, to check for flaring, or any periodic 
variation, but the object is too faint for any such variabil- 
ity to be noted. Since no sources are seen within 10" of 
this object, we used a region of radius 6" to extract the 
spectrum. As a measure of the PSF at the source location, 
a radial region of ~ 1.5" would enclose 50% of the energy 
at 1.49 keV. 

CXU J013623.4+154459 is very faint, and the higher 
energy {E > 9 keV) spectrum is seen to be strongly dom- 
inated by the background - possibly due to high energy 
particles, rather than X-ray events. Because of the prob- 
lems associated with accurately modeling the high en- 
ergy component of the background spectra, we used back- 
ground subtracted data and restricted our analysis to the 
energy range 0.3 — 10.0 keV. The background subtracted 
count rate in this energy band is 6.5(±2.4) x 10 -4 s _1 . Fix- 
ing the column density at 0.49 x 10 21 cm 2 , (obtained from 
online hydrogen column density generator of High Energy 
Astrophysics Archive (HEASARC), in the direction of the 
SN) and using the x 2 -Gehrels statistic, we found that the 
source was well fitted by a power-law spectrum with a 
photon index 0.6^ ; 35 , implying a 0.3-10.0 keV flux of 
1.3 x 10 -14 ergs cm -2 s _1 . The \ 2 goodness-of-fit param- 
eter was Q = 0.899, where Q measures the probability of 
obtaining the observed (or larger) value of x 2 , if the data 
was well represented by the fitted parameters. For an ideal 
fit, Q = l. 
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2.2. XMM-Newton observations of SN 2002ap field: 

SN 2002ap was observed by the XMM instruments from 
Feb. 2.03 to Feb. 2.42 UT, 2002 (estimated 5 days after the 
explosion), for a duration of 37.4 ksec. Both EPIC-MOS 
and EPIC-PN observations were carried out in the 'Thinl' 
filter mode. The data was pipeline processed and cali- 
brated using the XMM-Science Analysis Software (SAS) 
version 5.2, and the latest available versions of the calibra- 
tion files. The events in the EPIC-MOS 1, EPIC-MOS2 
and EPIC-PN data sets were filtered using the SAS- 
xmmselect task, using the appropriate flags and event se- 
lection criteria to account for event pile-up. We followed 
the analysis procedure recommended in "Users Guide to 
the XMM-Newton Science Analysis System" (2001)), with 
the exception of allowing patterns < 4 in the case of EPIC- 
PN, to account for any pile-up, however negligible it may 
be. A time series analysis of the entire EPIC-PN, MOS1 
and MOS2 data, shows that there are times of large, ran- 
dom fluctuation over the entire field of view of both EPIC- 
PN and MOS CCDs. Filtering out these intervals of highly 
fluctuating background reduced the exposure time to 25.5 
ks in EPIC-PN and - 30 ks in each of EPIC-MOS1 and 
MOS2. 

SN 2002ap was identified on the EPIC-CCDs using 
the optical coordinates (RA = 0l' l 36 m 23"85, DEC = 
+ 15°45'13"2). The EPIC-PN spectra was extracted us- 
ing a circle of radius 40" , within which almost 80% of the 
source photons in the energy range 10.5 — 1.5 keV would 
be enclosed, with the extraction region remaining on a sin- 
gle CCD chip. The EPIC-MOS1 and EPIC-MOS2 spectra 
were extracted using a psf of radius 50" , which encloses 
90% of the total energy. The response matrices and ancil- 
lary response files were generated using the SAS-rmfgen 
and SAS-arfgen tasks, and the appropriate line number for 
the PN-CCD. The data was grouped and analysed using 
the XSPEC software. 



2.3. Results of X-ray observations 

The 0.3-12 keV, background subtracted, count rate in 
XMM EPIC-PN is 6.6(±1.0) x 10~ 3 s" 1 . In the same en- 
ergy range, the EPIC-MOS1 count rate is 1.69(±0.64) x 
10~ 3 s- 1 and the EPIC-MOS2 count rate is 1.7(±0.64) x 
10~ 3 s _1 . 

In order to account for the presence of the Chandra 
source discussed above in the XMM spectral extraction 
region, all models fitted to the data had an additive ab- 
sorbed power-law component, with parameters derived 
from fits to the Chandra data, as discussed above. Since 
the SN is very faint, and the spectra rather sparse (Fig. 
[l]), we kept Nh — 0.049 x 10 22 cm -2 as constant for all 
fits. In general, we found that letting Ng vary can often 
result in unphysically low values of Nh <C 10~ 16 cm -2 , 
while providing little improvement in the value of reduced 
X 2 ■ In Table [l], we have quoted the parameters from the 



Table 1. Best fit spectral parameters for EPIC-PN. kT 
represents the plasma temperature or the black body tem- 
perature, depending on the model. All uncertainties are 
quoted in 90% confidence limit. Nh was constant for all 
fits. The absorbed flux for SN 2002ap is quoted here. 



Model 


N H 


a 


kT 


Xt/dof 


/ 0.3-10 




W 21 








10" 14 erg 




cm -2 




keV 




/cm 2 /s 


Power- law 


0.49 


Z - D 0.5 




1.2/20 


1.07 


Thermal 


0.42 
0.49 


o c0.6 
^•00.5 


0.84l° o j 


1.2/20 
1.2/20 


1.0 
0.81 


Brems. 












Raymond- 


0.49 




2-311JJ 


1.58/20 


1.04 


Smith 












Blackbody 


0.49 




0.2lt° o i 


1.4/20 


0.6 



best fitted models in our analysis - the 0.3-10 keV flux 
quoted there is corrected for the presence of the Chan- 
dra source. The data is rather sparse, hence x 2 -fitthig 
is unable to distinguish between thermal bremsstrahlung 
and the simple power law models. The power law dis- 
tribution (Fig. []]) fits the data well with a = 2.6, as 
does the thermal bremsstrahlung model with tempera- 
ture kT ~ 0.8 keV. Using the powerlaw model, we find 
that the total flux in the EPIC-PN extraction region is 
2.4 x 10~ 14 ergs cm" 2 s" 1 . Correcting for the presence of 
CXU J013623.4+154459, the 0.3-10.0 keV flux from from 
SN 2002ap is l.Q7±g;|f x 10" 14 ergs cm" 2 s" 1 . 

Adding an extra power-law, or a cutoff power-law 
component to the bremsstrahlung, or Raymond-Smith or 
blackbody models did not improve the fit, and resulted in 
poorly constrained values of temperatures with an unusu- 
ally high photon index a ~ 9.0, and hence these models 
were discarded. 

Because of the low count rate in EPIC-MOS CCDs, 
it was decided to fit to the combined MOS1 and MOS2 
data, allowing only the relative normalisations to vary. As 
discussed for EPIC-PN data, each model incorporates an 
absorbed power-law component to account for the Chan- 
dra source. Thus, the data was well fitted by an absorbed 
thermal bremsstrahlung, with x 2 /d.o./. = 0.95/19 and 
kT = 0.4l:^ 8 22 . The 0.3-10 keV flux (corrected for CXU 
J013623.4+154459 ) was 7.1 x 10" 15 ergs cm" 2 s _1 . 

SN 2002ap was too faint to be detected by either of 
the RGS detectors. 



2.4- The Optical Monitor observation. 

The XMM-Newton Optical Monitor observed SN 2002ap 
in the 4-frame (ENG-2) mode with UVW1 filter for a du- 
ration of 2.5 ks. The UVW1 filter peaks at A ~ 270 nm, 
and spans the range ~ 250 < A <~ 310 nm. The frames 
were combined and the data was pipeline processed using 
SAS version 5.3 with the latest available calibration files. 
A source list was compiled. We find that the optical mon- 
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pn_specgrpsrc.ds 




Table 2. 

brators 



GMRT observations for Flux and Phase Cali- 



channel enerav fkeVl 



Fig. 1. EPIC-PN spectrum and the fitted power law 
model (see Sect. |J) of SN 2002ap. The model fitted here 
consists of an absorbed powerlaw of photon index = -2.6 
with an additional constant absorbed power law compo- 
nent to account for CXU J013623.4+154459. \ogN H was 
fixed at 20.68. The residuals are plotted in the lower panel. 



itor count rate for SN 2002ap was 14.67(±0.08) s~\ with 



the background count rate being 0.8 s . This corresponds 
to a UVW1 flux of 7.667(±0.002) x 10 

A" 1 . 



-15 



erg cm 



3. GMRT observations of SN 2002ap field 

The radio observations of SN 2002ap were made with the 
Giant Metrewave Radio Telescope in aperture synthesis 
mode on 2002 February 5.56 in the 610 MHz band and on 
2002 April 8.22 in the 1420 MHz band. 

For both the observations at 610 as well as 1420 MHz, 
the pointing centre was chosen to be the position of super- 
nova i.e. RA= 01 ?l 36 m 23!85 and Dec= 15deg45T3'/2. In 
610 MHz observations the time on source were around 2 
hours and 3C48 was used as flux calibrator as well as phase 
calibrator. It was also used as bandpass calibrator. 3C48 
was observed for approximately 10 minutes after every 30 
minutes. 3C48 flux at 610 MHz was derived from the best 
fit spectra given in VLA calibrator manual. Out of a total 
of 30 antennas, we had 22 good antennas. The full band 
of 16 MHz was used, resulting in 128 frequency channels 
of width ~ 125 kHz. For the observations in 1420 MHz 
band, the number of antennas used were 24. The time on 
the source was ~ 3.5 hours. 3C48 was used as a flux cali- 
brator and it was observed at the beginning and at the end 
of the observations. The flux of 3C48 was assumed to be 
16.15 Jy from VLA calibrator manual. Compact sources 
0119+321 and 0238+166 were used as the phase calibrator 
with fluxes 2.6 Jy and 1.26 Jy respectively. 3C48 was also 
used as a bandpass calibrator. The full band of 16MHz 
resulting in 128 channels of width ~ 125 kHz was used. 



Source 


Date in 


RA 


Dec. 


V 


Flux 




2002 


J2000 


J2000 


MHz 


Jy 


3C48 


5 Feb 


01 37 41.3 


33 09 35 


610 


29.4 


3C48 


8 Apr 


01 37 41.3 


33 09 35 


1420 


16.2 


0119+321 


8 Apr 


01 19 35.0 


32 10 50 


1420 


2.9 


0238+166 


8 Apr 


02 38 38.9 


16 36 59 


1420 


1.0 



Table 3. Results of GMRT observations of SN 2002ap 



Date in 


u 


Resolution 


2a Flux 


RMS 


2002 


(MHz) 


(arcsec) 


(mJy) 


mJy/ beam 


5Feb 


610 


9.5 x 6 


< 0.34 


0.17 


8Apr 


1420 


8x3 


< 0.18 


0.09 



3.1. Data reduction, calibration and results 

The data was analysed using AIPS (Astronomical Image 
Processing Software) of National Radio Astronomy Ob- 
servatory. Interference in some of the frequency channels, 
spikes and dropouts due to electronics etc. were flagged 
using the AIPS task FLGIT. FLGIT was run on the cali- 
brator scans and 10 times the expected noise per channel 
was used as the flagging criterion. In this manner 15% of 
the data was flagged. Then one clean channel was looked 
at and antenna and baseline based flagging was done. Ap- 
proximately 20% data was lost in the flagging process. 
After computing the single channel antenna-based com- 
plex gains for the flux calibrator and phase calibrator with 
the time resolution of 20 seconds, bandpass calibration 
was done using 3C48 to apply the single channel calibra- 
tion to all the channels. To get rid of the effect of the 
band smearing at low frequency (610 MHz), a pseudo- 
continuum data base of 6 frequency channels was made 
from the central 100 channels of 16 MHz bandwidth. Field 
maps were made by Fourier transformation and cleaning. 
Iterative self-calibration rounds were run for phase correc- 
tion. Details of GMRT data analysis procedures are given 
in jBhatnagar, 200C| ). 

In neither of the observations carried out on 2002 
February 5 (day 8.56) and on 2002 April 8 (day 42), SN 
2002ap was detected by GMRT. The upper limits of the 
fluxes from the region of the SN are given in Table 3. 
The GMRT 610 MHz contour map shows several extended 
sources, which are likely associated with the spiral arms of 
the galaxy. Superposition of the 610 MHz contour map on 
the X-ray EPIC-PN image (Fig. ||), reveals the existence of 
4 X-ray sources with radio counterparts. These are listed 
in Table |4| with their radio flux densities and X-ray 0.3- 
10.0 keV count rates. Spectral analysis to determine the 
nature of these sources is being done separately. 
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Fig. 2. GMRT image of the field at 610 MHz containing 
SN 2002ap overlayed on the X-Ray image from XMM. 
X-ray image is on a grey scale while the radio image is 
the contour map with lowest level of 1.2 mJy. The image 
resolution is 15" x 15" . Four X-ray sources having radio 
counterparts are seen in this figure. They are numbered 
according to their decreasing RA in Table |[ Position of 
the SN on the map is marked with a cross. 

Table 4. List of X-ray sources with radio counterparts in 
M74. 



Source 


RA 


Dec. 


Radio 


EPIC-PN 


No. 






Flux 


count-rate 




J2000 


J2000 


(mJy) 


10" 3 ct s _1 


1 


01 36 47.2 


15 47 45 


7.1(±1.3) 


4.3(±1.3) 


2 


01 36 46.1 


15 41 17 


12.8(±1.5) 


3.9(±1.3) 


3 


01 36 24.9 


15 48 58 


22.7(±1.5) 


1.9(±1.3) 


4 


01 36 30.5 


15 45 17 


4.5(±1.3) 


4.6(±1.3) 



4. The location of the radiosphere 

The radio turn-on of supernovae are found to be 
wavelength dependent, with the flux at shorter wave- 
lengths peaking before that at longer wavelengths (e.g. 
Bergcr et al., 200"2| ) . This decreasing absorption can be ei- 
ther via free-free absorption (FFA) or synchrotron self- 
absorption (SSA) in the expanding CSM. We have fitted 
the VLA and GMRT data on 2002 Feb 5.96 to a SSA 
model, with spectral index = -0.8 in the optically thin 
limit, implying that the radius of the radio photosphere 



Table 5. Optical photosphere of SN 2002ap on day 5 
(Feb 2, 2002). The distance modulus was /x = 29.5, and 
the maximum value of bolometric correction was used 
( [Allen, 1976| ). 



M v 


M bol 


T eff 


Ropt 


Vph 


FjjVOIR 






K 


cm. 


km/s 


ergcm _2 s _1 


-17.4 


-16.5 


11000 


3.4 x 10 14 


8000 


2 x 10~ lu 



Table 6. Best fit parameters of SSA for SN 2002ap on 
day 8.96 i.e. Feb 5.96, 2002 



a Up 


F p 


Hr 


B 


GHz 




cm. 


G 


0.8 2.45 


397 


3.5 x 10 ib 


0.29 



on this day was R r — 3.5 x 10 15 cm and the magnetic field 
in the shocked ejecta was B = 0.29(5 (see Table ^|). We 
also note that that SSA prediction of flux at 610 MHz 
band is consistent with the GMRT upper limit. 

In order to compare the photosphcric radii across 
wavebands, we also compute the size of the optical pho- 
tosphere at the same epoch. Using the synthetic light 
curve best-fitted to the bolometric light curve on day 
5 QMazzali et al., 2002; ), 
Mboi 



the bolometric magnitude was 
-16.5, and the visual magni tude was M v = —17 .1 



(using a distance modulus of 29.5 ( Sharina et al.,1996|) ). 
This implies that the radius of the optical photosphere on 
day 5 was R opt = 3.4 x 10 14 cm, and corresponding to- 
tal flux across UV, optical and IR bands was Fijvoir = 
2 x 10 _10 erg cm~ 2 s _1 . The inferred conditions in the 
optical photosphere are listed in Table |. 

In contrast, the position of the radio photosphere on 
day 7 was r = 4.5 x 10 15 cm and its velocity 90, 000 km s _1 
( Berger et al., 200*2] ). Scaling to day 5, i? ra dio ~ 3.2 x 
10 15 cm. Thus, the radio photosphere is quite far outside 
the optical photosphere. The region of production of the 
X-rays is also outside the optical photosphere radius (see 
Sect. p). This substantiates the inference from the early 



spectrum (Meikle et al., 2002) that there was sufficient 
material at v > 30, 000 km s _1 . 

The framework of spherical symmetry has been used 
to perform the above analysis. This is a reasonable as- 
sumption on the day of the XMM observation and shortly 
afterwards, since the month-long, early spectropolarimet- 



ric data of Wang et al., 2002 implied there was little dis 



tortion of the photosphere on 2002 February 3. 

Finally, we compare the early radiospheric proper- 
ties of SN 2002ap with those of other radio SNe, across 
type classification, using equipartition arguments. The- 
oretically, early radio emission from SNe is thought to 
originate from an envelope that is expanding at a sub- 
s tantially highe r speed than the optical photosphere 
( 5hklovsky,1985 ). It is only at early times after the ex- 
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Comparison of spectrum of SN 2002ap, 1 993J and 1998bw on day 




fO s 10 1 " 
Frequency (Hz) 



Fig. 3. Comparison o f spectrum o f three SNe (SN1993J 
( |van Dyk et al., 1994|) , SN1998bw ([Kulkarni et al., 1 



and SN 2002ap ( |Berger et al., 2002| ) ) near day 11 after 
explosion. Solid line, dashed line and long dashed lines 
show the best fit synchrotron self absorption model for 
SN 2002ap, 1993J and 1998bw respectively and dotted 
line, dot-dashed line and short-long dashed line show the 
corresponding free free absorption model. 

Table 7. SN 1993J, SN 1998bw and SN 2002ap best fit 
parameters of SSA on day 11 



SNe 


v v 


F P 


e e , 


u eq 


Bo 


Ro 










xlO 45 




xlO 15 




GHz 


mjy 


fias 


erg 


G 


cm 


2002ap 


2.45 


0.48 


39.0 


0.69 


0.47 


4.80 


1998bw 


5.5 


50.4 


112.4 


3500 


0.23 


68.4 


1993J 


30.5 


22.3 


17.6 


0.50 


3.54 


1.08 



plosion that the effect of Compton boosting by electrons 
to high energy bands producing nonthermal spectrum is 
potentially observable (Sect. ^). At such early times radio 
spectra of supernovae are relatively rare. 

We have fitted models (Fig. ||) to the radio spectra ob- 
tained from three supernovae: SN 1993J (a type lib SN) 
on day 11.45, SN 1998bw (a type Ic SN) on day 11.7 and 
SN 2002ap (another type Ic SN) on day 11.0. Best fit pa- 
rameter values for SSA model were obtained by using a x 2 
fit. While SN 1993J and SN 1998bw were much brighter at 
their peaks compared to SN 2002ap, the synchrotron self- 
absorption peak for SN 2002ap seems to occur at the low- 
est frequency among all three. Best fit parameters values 
are listed in Table m. The equipartition (angular) radius 
(that gives the minimum energy density in combined mag- 



netic and rclativistic particle energy densities) obtained by 
|Slysh, 1963| : 



cq - 120d Mpc £> pmJy v PiGHz Mas 

appears to be the largest for SN 1998bw (112 /xas) while 
those for SN 1993J and SN 2002ap are comparable (18 and 
39 /ias respectively). Hence SN 1998bw which was not only 
coincident with a GRB source but also very bright in the 
radio was rapidly expanding. It was unusual compared to 
the more ordinary SN like SN 2002ap (despite both su- 
pernova's identical type classification), - the latter hav- 
ing similarity with SN1993J in terms of equipartition size, 
across type classification. 



5. Radio &: X-ray turn-on & progenitor mass loss 

The early emission in X-ray and Radio bands allows us to 
estimate the density of the CSM along line of sight, and 
hence the mass loss rate from the presupernova star, if 
the terminal wind speed of the progenitor star is known 
or estimated. 

The CSM density is related to the mass loss param- 
eters as: p cs — M /A-Ku w r 2 . The gas overlying the emis- 
sion regions would itself serve as an absorber of the ra- 
dio and X-ray radiation. An estimate of M can be made 
from the XMM data. The mass loss rate is a measure of 
the column depth J ndr rather than the J n 2 dr which ap- 
pears in the radio (free-free) absorption. The stellar wind 
would absorb the X-rays below 8 keV due to the dom- 
inant photo-electric absorption of the metals. The time 
at which t *~ 1 is reached for photons of energy -EkoV is 
(|Chevalier fc Fransson, 1994[): 



<x = C 5 M_ 5 u~ 



,-lp-8/3 
4 -^keV 



The value of the constant C5 (n) depends upon the den- 
sity profile (n) of the ejecta (p C j = p (t/t )~ 3 (ut/r) n ). 
Using a typical value of the profile index n = 10.2, the 
epoch of X-ray detection as day 5 and Ek e v = 0.7, we 
determine A/_ 5 /v wl < 6.5 x 10~ 2 . From optical spectral 
measurements, the velocity of the pre-supernova wind was 
~ 580 km s" 1 ( |Leonard et al., 2002| ). This implies that 
M < 4 x 1O _5 M . The above expression assumes that 
the heavy elements contributing to the opacity are mostly 
neutral - this upper limit would be less stringent if Oxy- 
gen in the undisturbed circumstellar medium is heavily 
ionised, in which case the X-ray absorption is reduced. 

A similar mass loss rate can also be derived from the 
first radio detection of SN 2002ap in the 1.4 GHz band 
on 2002 Feb 1.93 (i.e. on day 5) when the optical depth 
in the radio frequency band would have been ~ 1. At this 
time the dominant opacity source for radio radiation is the 



free- free opacity from fully ionised wind (Chevalier, 1982): 
From the radio data, we estimate the upper limit on the 
mass loss rate as M < 6 X 10 _5 Mq yr" 1 . 
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6. Radiative processes in the early SN 

Our spectral modeling of XMM observation of SN 2002ap 
shows that both power-law and thermal bremsstrahlvmg 
models give acceptable fits to early X-ray data. Since mul- 
tiwaveband data from UV, optical and radio are also avail- 
able from this supernova, these can be additionally utilised 
to restrict models of prompt emission. Among the possi- 
ble mechanisms invoked in the context of radio loud X-ray 
supernovae are: a) direct synchrotron emission from radio 
to X-ray band, b) free-free emission from the hot gas con- 
sisting of swept up mass behind the blast wave (circum- 
stellar) shock and that behind the reverse shock, and c) 
inverse Compton scattering of optical photons emitted by 
the supernova by the hot electrons. 

Case (a): As in the case of SN 1980K 
(Canizares et al, 1982), the direct synchrotron radi- 
ation mechanism for SN 2002ap prompt emission 
can be eliminated. The power law spectrum in fre- 
quency for the radio flux observed from SN 2002ap is: 
F v = 3mJy(3.2GHz/^) a9 (t/day)-°' 9 (unabsorbed flux 
from Berger et al., 2002). If the same ambient magnetic 
field (B ~ B eq ~ 0.2 Gauss) implied by equipartition 
arguments were to hold, to generate X-ray energies 
(1 keV) by synchrotron mechanism, would imply the 
presence of highly relativistic electrons (7 > 1.4 x 10 6 ). 
On day 5, a straight extrapolation of the radio spectrum 
with the above flux would imply a flux of only 58 
picoJansky at 1 keV. Since the 90% energy bandwidth 
for a power law distribution of energy index -1.6 in the 
XMM band is: 3.7 keV, this would imply an X-ray flux 
of about 5 x 10~ 16 erg s _1 cm -2 . This is far less than 
what is observed by XMM. Thus a direct radio to X-ray 
synchrotron radiation mechanism is unviable, or produces 
an insignificant part of the observed X-rays. Therefore, 
different populations of electrons appear to be responsible 
for the radio and X-ray fluxes that are seen. 

Case (b): The circumstellar gas is heated to high 
temperatures (> 100 keV) by the shock resulting 
from the piston-like expanding supernova envelope. This 
shocked gas cools by free- free emission and/or Comp- 
ton cooling. The free- free luminosity from the circum- 
stellar (blast- wave) and reverse shocks (Fransson, 1982, 



Chevalier & Fransson, 2001) is: 



U = 3 x 10 jy 5ff C n 



M_5 



t 



10 day 



erg s 



where for the reverse shock, C n = (n— 3)(n— 4) 2 /4(n— 2) = 
7.875 for n=10 and C n — 1 for the circumstellar shock, gff 
is the Gaunt factor. 

With M = 3.36 x 10 _5 M Q yr" 1 and with a wind ve- 
locity u w — 1000 km s _1 (i.e. u w i = 100), at t=5 days, 
the X-ray flux at 7.3 Mpc would be: 8.2 x 10~ 15 erg cm -2 
s _1 . This flux would have been near the observed flux 
that has been detected from SN 2002ap by XMM. How- 
ever, with the reported high ejecta velocities and the im- 



plied temperature of the shocked ejecta and CSM along 
with the limited cool absorbing shell at this early stage, 
one would expect a flat tail of high energy photons up 



to about 100 keV ( [Fransson, 1982[ ). Since the XMM spec- 
trum on day 5 is not so hard (see Table 1 where the ther- 
mal bremsstrahlung temperature is quite modest), the free 
free emission is also an insignificant part of the early radia- 
tive budget. Thermal X-ray emission (i.e. bremsstrahlung) 
will however become more dominant compared to the in- 
verse Compton process (case c below) as the supernova 
ages. 

Case c): The dominant cooling and radiative mech- 
anism at early stages when the optical photospheric 
temperature T e f f > 10 4 K, is Compton cooling 



( Fransson, 1982 ). Optical photons from the photosphere 
can undergo repeated scattering off the hot electrons 
in the shocked region (notably the shocked CSM) and 
a power law photon distribution (in energy) can result 
even if the electron distribution is not a power law, and 
even if the electron scattering depth r c remains small. 



Pozdnyakov et al., 1977 reported Monte Carlo calcula- 



tions and their analytical approximations of the emergent 
Compton scattered spectrum. The spectrum depends on 
the optical depth r c and the temperature T of the electron 
plasma. For photons undergoing electron scattering with 
cross-section <jt in the shocked gas is: 



Mer T 



( 1 _ ^opt \ 

v R s ) 



47TTO D i? s U, 



The flux density and the energy index 7 originating from 
Comptonization are ( jFransson, 1982 ) : 



with 



7(7 + 3) 



m c c 



In 



-(0.9228 -lnr e ) 



Here, v a is the peak frequency in the (optical) input spec- 
trum, and F Q the optical flux density. i? op t is the radius 
of photosphere. R s is the position of the plasma element. 

We note that on day 5 the unabsorbed (or dereddened) 
X-ray and optical flux densities derived from XMM and 
ground based observations (Sect. ^) imply an optical to 
X-ray power law of 7 = 2.5 — 2.8 and a logarithm of the 
ratio of flux densities of « 7.4. The optical to X-ray power 
law index is consistent with the XMM results (see Table 
0) . We determine the optical depth and temperature con- 
ditions of the Comptonizing plasma for 7 = 3 in Table 
8. The presupernova scenarios typical for type Ic SNe, for 
which the optical depths are reported in Table || are dis- 
cussed in the next section. The plasma has the maximum 
optical depth at twice the optical photosphere radius. The 
latter was taken as 3.4 x 10 14 cm. Since most of the X-ray 
emission would take place at ~ r max , the relevant plasma 
outflows with a velocity of approximately 16, 000 km s . 
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Table 8. Comptonizing plasma properties at t = 5d with 
7 = 3 and R opt = 3.4 x 10 14 cm 



Scenario 


M_ 5 


U w i 




T e 




lO- 5 M /yr 


10 km/s 




10 9 K 


Wolf-Rayet 


1.5 


58 


4 x 10" 4 


2 




3 


100 


4.4 x 10" 4 


2 


Interacting 


10 


58 


2.5 x lCT 6 


1.5 


Binary 










Case-BB 


10 


10 


1.5 x 10" 2 


1.1 



7. Discussion and Conclusion 

In this paper, we presented an analysis of the XMM data 
of the SN 2002ap field and have obtained spectral model 
fits to the prompt X-ray emission. We compare the X- 
ray image with the GMRT 610 MHz radio image obtained 
three days apart. While we find no radio counterpart of 
the SN at such low frequencies, several sources in the field 
have radio and X-ray counterparts. 

We compare the radio data obtained from three differ- 
ent supernovae in their early phases and model these us- 
ing the synchrotron self absorption model. SN 1998bw(Ic) 
with a GRB counterpart had very different radiosphere 
radius and equipartition angular sizes at approximately 
the same time in their evolution compared to two other 
SNe: SN 2002ap(Ic) and SN 1993J(IIb). 

We model the early X-ray emission with inputs from 
optical photometry and light curve and find that the in- 
verse Compton scattering of optical photons from the su- 
pernova photosphere by hot electron plasma can account 
for the observed early X-ray flux and the spectrum for 
modest electron temperatures and optical depths. Ther- 
mal processes are inefficient initially and would be impor- 
tant only as the supernova ages. 

Mass loss rates and stellar wind velocities of the pro- 
genitor stars determine the optical depth of shock heated 
matter due to electron scattering, - a key parameter in 
the production of the X-ray flux from the lower energy 
photons. These in turn depend upon the scenario and pro- 
genitor configuration giving rise to type Ic SNe, e.g.: 

1. Massive Wolf-Rayet (WR) stars which have lost 
their hydrogen and helium envelope before the explosion 
( Langer et al., 1999 ) have an empirically determined mass 
loss rate of 1.5 — 3.2 x 10~ 5 M Q (afte r taking into account 
effects due to clumpy medium, ( |Hamann et al., 1998 , 



Willis, 1998 ) and terminal wind velocities ~ 1000 km s" 
depending upon the type of the WR star. 

2. Interacting binaries - in particular, in a case BB 
mass transfer from a helium star overflowing its Roche- 
lobe to a companion removes most of the helium rich lay- 
ers before the type Ic SN. Habets (1985) finds that a 2.5 
Mq helium star during carbon shell burning stage expands 
to a red giant dimension of 18 Rq and sustains an aver- 
age mass transfer rate of 10 _4 MQyr _1 lasting about 3000 



years. During this time, the terminal wind speed, of the 
mass losing star would be typically 100 km s -1 . 

Using the relevant parameters in the above two sce- 
narios, we derive electron optical depth encountered by 
the intermediate energy photons as listed in the Table |[ 
It is evident that both ranges of electron optical depths 
for the Compton boosting process remain viable alterna- 
tives. For the case of the interacting binary model where 
the optical depths are somewhat larger, the implied elec- 
tron temperatures required for the plasma would be lower 
than in the single star WR model. Such temperatures 
are well within the range expected for hot circumstel- 
lar gas (T cs > 3.6 x 1Q 9 K) even for modest velocities 
of 16,000 — 20,000 km s _1 for the hot electron plasma 
moving above the optical photosphere. 
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